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Abstract 

Production of even-parity charmed mesons in hadronic B decays is studied. Specifically, we focus 
on the Cabibbo-allowed decays B — ► D**ir and £)**£)(*), where D** denotes generically a p-wave 
charmed meson. While the measured color-allowed decays B° — > D** + tt~ are consistent with the 
theoretical expectation, the experimental observation of B~ — > D**°7r~ for the broad D** states is 
astonishing as it requires that the color-suppressed contribution dominates over the color-allowed 
one, even though the former is l/m& suppressed in the heavy quark limit. In order to accommodate 
the data of B — ► D**tt~ , it is found that the real part of 02/01 has a sign opposite to that inB-> Dn 
decays, where oq and ai are the effective parameters for color-allowed and color-suppressed decay 
amplitudes, respectively. The decay constants and form factors for D** and the Isgur-Wise functions 
Ti/ 2 (uj) and t 3 / 2 (w) are extracted from the data of B — > Z)**7r decays. The Isgur-Wise functions 
calculated in the covariant light-front quark model are in good agreement with experiment. The 
neutral modes B° -» D**°tt° for D** = D$ (2400) , D[ (2430) and B° -> Dj° (2430)cj are predicted 
to have branching ratios of order 10 -4 which are also supported by isospin argument. The decay 
constants of D* (2317) and D' sl (2A6Q) are inferred from the measurements of B — > D**~D to be 
58 ~ 86 MeV and 130 ~ 200 MeV, respectively. Contrary to the decay constants fp* and fpi 
which are similar in size, the large disparity between fjj* and fjj> is surprising and unexpected. 
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I. INTRODUCTION 



The spectroscopy for the p-wave charmed mesons has entered a new and exciting era since 2003. 
First of all, BaBar |l| has discovered a new narrow and light resonance -D* (2317). The existence of 
a second narrow resonance Z) s j(2460) which can be identified with J p = 1 + state was first hinted 
by BaBar Q| and then observed and established by CLEO 2] and Belle 3]. Second, the broad 
Dq and D\ resonances which are the counterpart of D* (2317) and -0^(2460) in the non-strange 
charm sector, were explored by Belle Q] in charged B to D + 7r~7r~ and D* + tt~it~ decays and by 
FOCUS 5<] in photoproduction experiment. 

Just like the light scalar mesons /o(980), ao(980), c(600), the underlying structure of the above- 
mentioned p-wave charmed mesons is not well established theoretically. Recall that before 2003, the 
p-wave states D* and D' sl with j q = 1/2 (j q being the angular momentum of the light degrees of 
freedom) are predicted to be broad and decay into DK and D*K, respectively, in the conventional 
quark model [6(. However, the observed -D* (2317) and -0^(2460) states are below the DK and 
DK* thresholds, respectively, and hence are very narrow. This unexpected and surprising disparity 
between theory and experiment has sparked a flurry of many theory studies. It has been advocated 
that this new state is a DK molecular [?J or a D s tt atom 0] or a four-quark bound state (first 
proposed in 0], followed by 0, On the contrary, it has been put forward that, based on 

heavy quark effective theory (HQET) and chiral perturbation theory, the newly observed -D s (2317) 
is a + cs state and that there is a 1 + chiral partner with the same mass splitting with respect to 
the 1~ state as that between the + and 0~ states 0, namely, m D >^ — rriD* = "t,d* — m D a - 

The spectra and strong, radiative decays of the p-wave charmed mesons have been studied 
extensively (for a review, see 

The measurements of two-body B decays into D** and 
D**, where D** denotes a generic even-parity charmed meson, provide very useful information on 
the decay constants and form factors of the excited charmed meson and the ratios of radiative to 
hadronic decay rates, for example, T(D* Q — > D*'j)/F(D* — » D s ir°). Moreover, they provide an 
opportunity to test heavy quark effective theory. As we shall see below, the decay amplitudes of 
B — > D**tt in the heavy quark limit are governed by the Isgur-Wise functions Ti/ 2 (^) and r 3 / 2 (w). 

In the present work, we will study even-parity charmed meson production in B decays. Specifi- 
cally, we focus on the Cabibbo-allowed decays B — > D**tt and D**D^*\ 2 The decay B — > D**°ir~ 
receives color-allowed and color-suppressed contributions, characterized by the effective Wilson pa- 
rameters a\ and (Z2, respectively. It is important to know whether the interference between the a\ 
and a2 terms is constructive or destructive. In the naive factorization approach, a\ = c\ + C2/N c 
and 02 = C2 + c\/N c . In the late 80's and early 90's, the large- N c approach was very popular and 



1 An issue for the four-quark model is whether there exist the conventional cs and cq states. A non- 
observation of a heavier and broad + cs state will not support the four-quark interpretation of Z?* (2317). 
For non-strange charmed mesons, it has been argued that Dq(2308)° observed by Belle Q is a four-quark 
state, whereas -Dg(2405)°' + measured by FOCUS is a normal eg state that is, the Dq state observed 
by Belle and FOCUS may not be the same. 

2 These decays have been studied previously in[l3[ll[iaH|2ll|2jl0 
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successful in explaining the hadronic weak decays of charmed mesons 28(]. It was widely believed 
by most practitioners in the field that the 1/N C expansion applies equally well to the weak decays 
of B mesons. Since 02/(11 ~ C2/C1 ~ —0.25 at the renormalization scale fi = tub in the leading 
1/N C expansion, P° — > D + ix~ is naively expected to have a larger rate than B~ — » D°ir~ due to 
the destructive interference in the latter. However, the CLEO measurements of B — > Dtt imply 
the opposite 29] ! This is a very stunning result. In order to accommodate the B — > Dtt data, 
the ratio 02/01 is found to be order of 0.20 ~ 0.25. In early 2000's, the color-suppressed mode 
B — > D°ir° is found to be significantly larger than the theoretical expectation based on naive 
factorization. For example, the measurement B{B^ — > D°Tt°) = (2.91 ± 0.28) x 10 -4 30] is larger 



than the theoretical prediction, (0.58 ~ 1.13) x 10 |3J, by a factor of 2 ~ 4. Moreover, the three 
B — > Dn amplitudes form a non-flat triangle, indicating nontrivial relative strong phases between 
them. As a consequence, 02/01 ~ (0.45 — 0.65)6^°° is inferred from the B —* Dtt measurements 
including the neutral mode B° — > D°ir° 0, 0, E3 ■ The question is then why the magnitude 
and phase of Qs/qi are so different from the model expectation. To resolve this difficulty, it has 
been shown in [351 136 1 that the enhancement of the effective 02 and its strong phase can be ascribed 



to final-state interactions. 

For B — > D**7T decays, we will pay a great attention to the relative sign of the decay constants 
and form factors of the p-wave mesons. It turns out that in order to explain the larger rate of 
B - D **o n - than D **+ n - f or D ** = D * and D >^ trie real part of o 2 /oi for B -» D**ir 

has to be negative with a large magnitude. Moreover, the color-suppressed contribution has to 
dominate over the color-allowed one, in contrast to the naive expectation that the color-suppressed 
amplitude is 1 /mj suppressed in the heavy quark limit. This is a third surprise for the ratio 02/01 
! The question to be addressed is why the sign of the real part of 02/01 flips when B — > Dtt is 
replaced by B — > D**n. 

This work is organized as follows. In Sec. II, the decay constants of D** and B — > D** form 
factors within the covariant light-front quark model are summarized. The decays B — > D**tt and 
B — > D* S *D are studied in Sees. Ill and IV, respectively. Conclusions are presented in Sec. V. 



II. DECAY CONSTANTS AND FORM FACTORS 

In the quark model, the even-parity mesons are conventionally classified according to the quan- 
tum numbers J,L,S: the scalar and tensor mesons correspond to 2S+1 Lj = 3 i"b an d 3 -P2 5 re- 
spectively, and there exit two different axial- vector meson states, namely, l P\ and 3 Pi which can 
undergo mixing if the two constituent quarks do not have the same masses. For heavy mesons, the 
heavy quark spin Sq decouples from the other degrees of freedom in the heavy quark limit, so that 
Sq and the total angular momentum of the light quark j are separately good quantum numbers. 
The total angular momentum J of the meson is given by J = j + Sq with S = s + Sq being the 

total spin angular momentum. Consequently, it is more natural to use 

1/2 

Pq to classify the first excited heavy meson states where L here is the orbital angular momentum 
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TABLE I: The 


masses 


and decay widths of even-parity charmed mesons (30I] . 


The four p-wave 


charmed meson 


states 


are denoted by Dq,D[, D\ and 


D|- in the heavy quark limit, .D^ has 


j = 1/2 and D\ has j 


= 3/2 with j being the total angular momentum of the light degrees of 


ii eeciom. 










State 




Mass (MeV) 




Width (MeV) 


Z?o(2400)° a 




2352 ± 50 




261 ± 50 


-DqI^oo)* a 




2403 ± 14 ± 35 




283 ± 24 ± 34 


Di(2420) 




2421.8 ± 0.8 




on o + l 9 

20-3Ii.s 


Di(2420) 




2427 ± 5 




26 ± 8 


1)1(2430)* 




2427 ± 26 ± 25 




3841 ± 74 


DX(2460)° c 




2462.7 ± 0.8 




44 ± 2 


,D|(2460) ± 




2464.9 ± 3.0 




29 ±5 


D s * (2317) 




2317.3 ± 0.4 ±0.6 


i 


< 3.8 d 


£^(2460) 




2458.9 ± 0.9 




< 3.5 d 


L> s i(2536) 




2535.35 ±0.31 




< 2.3 


D s * 2 (2573) 




2573.5 ± 1.7 




15±S 



"While the mass and the width of Dq(2400)° arise from the average of Belle |4| and FOCUS .5] measure- 
ments, the mass and the width of L»5(2400) ± are solely due to FOCUS [5j. 

including the most recent CDF measurements m(D?) = 2421.7±0.7±0.6 MeV and T(D%) = 20.0±1.7±1.3 
MeV 38] to the 2005 PDG average 0. 



c Including the most recent CDF measurements m(-D| ) 
2.3 ±1.3 MeV to the 2005 EDG average 0. 
d The width limit from BaBar 1391. 



2463.3 ± 0.6 ± 0.8 MeV and T(D 



49.2 ± 



of the light quark. It is obvious that the first and last of these states are 3 P 2 and 3 P , while |3J 



\p? /2 ) = ] fl\ i Pi)+ ] fl\ 3 Pi), \pi ,2 )=-sl\\ l Pi)+^m- (2.1) 



In the heavy quark limit, the physical eigenstates with J p = 1 + are P^ 2 and P^ 1 " rather than 6 P\ 
and 1 P 1 . 

The masses and decay widths of even-parity (or p-wave) charmed mesons Dj and D*j are 
summarized in Tabled It is known that in the non-charm scalar meson sector, the quark model 
cannot explain why the scalar strange meson .K"q(1430) with a mass 1412 ± 6 MeV [30j is lighter 
than the non-strange one ao(1450) with a mass 1474 ± 19 MeV 30]. 3 Likewise, it is clear from 
Table [I] that the relation m(D**) > m(D**) holds except for D* and Dq for which we have 
m(D* ) < m{Po). This is the place where the conventional quark model seems not to work 4 and 



1/2 



3 Recently it has been advocated in |4(j that the puzzle with the relative masses of Kq (1430) and ao(1450) 
can be resolved provided that the observed scalar nonet in the mass range 1-2 GeV is a tetraquark nonet 
plus a glueball. The yet-to-be observed qq scalar nonet lies around 1.1 GeV. 

4 There are some attempts in understanding the mass relation m(D* ) < iti(Dq) by modifying the conven- 
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a four-quark structure for Dq and D* is preferable 

We shall use l' + and 1 + or D[ and -Di to distinguish between two different physical axial-vector 
charmed meson states. 5 The physical 1 /+ state is primarily p\^ 2 , while 1 + is predominately P^ 2 . 
This is because, in the heavy quark limit, the physical mass eigenstates D[ and D\ can be identified 

1/2 3/2 

with P 1 1 and P l , respectively. However, beyond the heavy quark limit, there is a mixing between 
p\l 2 and P^ 2 , denoted by d\^ 2 and D^ 2 respectively, 

Di(2430) = ^ /2 cos6» + Di /2 sin6», 

£>i(2420) = -Dl /2 sine + Dl /2 cosO. (2.2) 
Likewise, for strange axial-vector charmed mesons 

D' sl (2460) = Dl{ 2 cos 6 a + D 3 s { 2 S mO a , 

Z? s i(2536) = -Dl{ 2 sin 9 S + D 3 S { 2 cos 9 S . (2.3) 

1/2 3/2 

Since -D/ is much broader than , the decay width of -Di(2420) is sensitive to the mixing angle 

1/2 3/2 

6>. The D{' —D] 1 mixing angle was reported to be 

= 0.10 ± 0.03 ± 0.02 ± 0.02 rad = (5.7 ± 2.4)° (2.4) 

by Belle through a detailed analysis of B — > D*ittt |4|. 

Since the decay B~ — > D**°ir~ receives color-allowed and color-suppressed contributions, it 
is important to know whether the interference is constructive or destructive. In [4ia] we have 
computed the decay constants and form factors for the ground-state s-wave and low-lying p-wave 
mesons within the framework of a covariant light-front approach. 6 In our approach, we first 
fix the vertex functions (i.e. Feynman rules for the meson-quark-antiquark vertices) for both s- 
wave and p-wave mesons. Then we are able to compute their decay constants and form factors. 
Hence, the relative sign and the factors of i between two-body and three-body matrix elements 
can be determined. We then adopt two different approaches to elaborate on the heavy quark limit 
behavior of physical quantities: one from top to bottom and the other from bottom to top. In the 
top-to-bottom approach, we derive the decay constants and form factors in the covariant light-front 
model within HQET and obtain model-independent heavy quark symmetry (HQS) relations. In 
the bottom-to-top approach, we study the heavy quark limit behavior of the decay constants and 
transition form factors of heavy mesons and show that they do match the covariant model results 
based on HQET |43j. 



tional potential model. For example, one loop chiral corrections to the potential model is considered in 
|4l|. while the potential model in 01 takes into account negative energy states of a heavy quark in a 
bound state. 

5 The notation for D[(2430) and L>i(2420) is opposite to that in |3l| . 



6 There are many typos in the printed version of |43| but not in the archive version, hep-ph/0310359 
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A. Decay constants 



The decay constants of scalar and pseudoscalar mesons are defined by 7 

(Q\Ap\P(q)) = ifpqp, <(TO5(g)> = feq^. (2.7) 

It is known that the decay constants of non-charm light scalar mesons are smaller than that of 
pseudoscalar mesons as they vanish in the SU(3) limit. The decay constants of the axial-vector 
charmed mesons are defined by 



( Q \ A n\ D l (?> e )> = / D i/2m D i/2^, (0\A^\Dl /2 (q,e)} = f ^3/2 m D s/2£ M . (2. 
It is known that in the heavy quark limit 44], 



/ n l/2=/ Dn *, / D 3/2=0. (2.9) 



Since the decay constant of D% vanishes irrespective of heavy quark symmetry (see below), the 
charmed mesons within the multiplet (0 + , 1 + ) or (l' + , 2 + ) thus have the same decay constant. 
The polarization tensor e^ u of a tensor meson satisfies the relations 

e^ = e^, 8^ = 0, = = 0. (2.10) 

Therefore, 

(0\(V - A)^\D*(e,p)) = oe^jT + be\ Vtl = 0. (2.11) 

The above relation in general follows from Lorentz covariance and parity considerations. Hence 
the decay constant of the tensor meson vanishes; that is, the tensor meson D\ cannot be produced 
from the V — A current. 

Using f D = 200 MeV, f Da = 230 MeV and f D * = 230 MeV as input, the decay constants (in 
units of MeV) of p-wave charmed mesons are found to be |43J 

Jd = 86, f D i/ 2 = 130, f D s/2 = -36, 

f Ds0 = 7l, / i /2 = 122, / 3/2 = -38, (2.12) 



7 Sometimes the decay constant of the pseudoscalar meson is defined as (0\Au\P(q)) = fpq^ in the literature 
(see e.g. This corresponds to choosing a non-hermitian vertex function r = 75 for the pseudoscalar 

meson or redefining the phase of the pseudoscalar field, namely \P) — ► exp(i<ft)\P) with <p = n. The 
two-body matrix element (P\V^\B) is still given by Eq. (|2.15|l . However, the usual soft-pion theorem 



.V2 



Hm<7r(«)|^|B) = -t^-(0\A^\B), (2.5) 



has to be modified to 



V2 



lim<7r(«/)|VgB) = ^(0\A^\B). (2.6) 
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in the covariant light-front model, where we have used the constituent quark masses 

m M = 0.26GeV, m s = 0.37GeV, m c = 1.40GeV, m b = 4.64GeV. (2.13) 

Notice that although D s has a decay constant larger than that of D as expected, it is other way 
around for the scalar mesons, namely, /d > /d* - This can be seen from the light-front quark 
model expression j^] 

/d s (d* ) oc / dx 2 ■ ■ ■ [m c x 2 ± m s (l - x 2 )\. (2.14) 

Since the momentum fraction x 2 of the strange quark in the D S (D* ) meson is small, its effect being 
constructive in the case of D s and destructive in D* is sizable and explains why fjj* //d s ~ 0.3 
and f D * o < f Do . 

In principle, the decay constants of the p-wave strange charmed meson D** can be extracted 
from the hadronic decays B — > DD** since they proceed dominantly via external VF-emission. In 
Sec. IV we shall extract the decay constants of D* and D' sl from experiment. 

There are other model calculations of the p-wave charmed meson decay constants. In general, 
these estimates are larger than ours, (|2.12j) . For examp le, the QCD sum rule approach in |45| 



yields f D * = 170 ± 20 MeV, while the quark model in |4<j predicts f D * = 139 ± 30 MeV and 
/d* — 110 ± 18 MeV. As we shall see below, some of the decay constants for p-wave charmed 
mesons can be phenomenologically extracted from B —* D**ir and B — > D**D decays and compared 
with model predictions. It turns out that while our prediction for fp* is smaller than experiment, 
our result for fu* is in agreement with the data. 



B. Form factors 

Form factors for B — » M transitions with M being a parity-odd meson are given by 
(P( P )\V,\B( PB )) = ((PB+^-^^S^ 

ly{ Pl e)\V,\B{ PB )) = \ e^ a pe*»p%pPv BV (q 2 ), 

m B +my 

{V(p,e)\A,\B( PB )) = i{(m B + m v )e;A? v (q 2 )- E * '* B (p B + p)^ '(q 2 ) 

<~ p m B + my 

-2m v £ -f- g^fV) ~ ^V)]}, (2-15) 

where e i23 = 1, q = p B - p, F BP (0) = F BP (0), A BV (0) = A BV (0), and 

ABV(q2) = rn^y AfV{q2) _ rr^ny ^ (?2) ^ 



For B — > D** transitions, we use 



(D* (p)\An\B(p B )) 



'(PB+P),- " 1 2 ml ° q ,)F BDo ( q 2 ) 



m B m D w BD Q(n 2\ 
— 2 fyA) (? ) 
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(D 1 (p,e)\V lt \B(p B )) = -iU mB -m Dl )elV^{q 2 )- ' PB ( P b + p)„V 2 BD ^q 2 ) 

{ ^ m B - m Dl 

£ ■ PB , . r fiflj, 2\ TrBDi ( 1 

- 2m Di q2 (PB ~ P)p [V 3 (q ) - V (q 

(D 1 (p,e)\A fl \B(p B )) = W£ ^/i BDl ( g 2 ), (2.17) 

m B — m/jj 

{D* 2 {pe)\V^\B(p B )) = h{q 2 )e^ pa e wa {p B ) a {p B +p)P{p B -pY, 
(D* 2 (p,e)\A^\B(p B )) = -i[k{q 2 )e%p B + b + (q 2 )e* ap p%p B {p B +p)^ 

+ b-(q 2 )e* al3 p B p B (p B - p)^ 



with 



lm B . l ^rno 1 



and Vf Dl (0) = V BDl (0). 

Note that except for the form factors h, all the other form factors are dimensionless. In 

principle, it is better to parametrize the form factors in such a way that they are all positive-defined. 
This is the case for B to s-wave meson transitions, but not for all B — > D** transition form factors. 
At any rate, the signs of various form factors can be checked via heavy quark symmetry shown 
below in Eq. (|2.21|) . For example, a factor of (— i) is needed in B — > S transition in order for the 
B — ► Dq form factors F^q° to be positive. 

Given the Feynman rules for the meson-quark-antiquark vertices (see Table I of in ^ ne 
framework of the covariant light-front (CLF) quark model, we are able to compute the form factors 
in the spacelike momentum transfer q 2 < 0. Form factors at q 2 > can be obtained by first 
recasting them as explicit functions of q 2 in the spacelike region and then analytically continue 

BD 3/2 

them to the timelike region. We find that except for the form factor V 2 , the momentum 
dependence of form factors in the spacelike region can be well parameterized and reproduced in the 
three-parameter form: 

F(?2) = l-a(qym 2 B )+b(qym 2 B y ' (2-19) 

for B — > M transitions. The form factor V 2 1 approaches to zero at very large —\q \ where the 
three-parameter parametrization (|2.19|) becomes questionable. To overcome this difficulty, we will 
fit this form factor to the form 

F{q2) = (1 - q 2 /m 2 B )[l - a(qym 2 B ) + 6((7 2 /m|) 2 ] (2 ' 20) 

and achieve a substantial improvement. 

Form factors for B — > ir and B — ► D** transitions calculated in the CLF model are listed in 
Tables ITTI For comparison, form factors evaluated in the Isgur-Scora-Grinstein-Wise (ISGW) quark 
model |48j are also exhibited in Table IIIII Before our work, the ISGW quark model is the only 
model that can provide a systematical estimate of the transition of a ground-state s-wave meson to 
a low-lying p-wave meson. This model is based on the non-relativistic constituent quark picture. In 
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TABLE II: Form factors for B — > tt, D** transitions obtained in the covariant light-front model [43j 
and fitted to the 3-parameter form Eq. (|2.19|) except for the form factor V 2 denoted by * for which 



the fit formula Eq. ([2.20)1 is used. All the form factors are dimensionless except for h, 6+, 6_ with 
dimensions GeV -2 . 



F 


F(0) 




a 


b 


F 


F(0) 




a 


b 


rpBw 
^1 


u.zo 


1 1 R 
1 . ID 


1 71 
1 . / 


u.yo 




u.zo 


U.oO 




n 1 n 
u. 1U 


F l 


n O/i 
0.Z4 


n A 
0.34 


1 no 
1.03 


n 0^7 
0.2 / 




n /I 
0.24 


n on 
0.20 


n /I n 
— 0.49 


0.35 


A BD i 2 


n 10 
— U. 1Z 


n 1 a 


n 71 

U. 1 1 


n 1 s 
u. 10 


BD \' 2 


U.Uo 


u. 10 


1 98 
l.Zo 


n 9Q 
— u.zy 


Rn l/2 

y BD x 


-0.19 


-0.13 


-1.25 


0.97 


bd\' 2 

v 2 


-0.12 


-0.14 


0.67 


0.20 


a bbT 


0.23 


0.33 


1.17 


0.39 


bdI /2 


0.47 


0.70 


1.17 


0.03 


v BD-r 


0.55 


0.51 


-0.19 


0.27 


bd\ /2 


-0.09* 


-0.17* 


2.14* 


4.21* 


h 


0.015 


0.024 


1.67 


1.20 


k 


0.79 


1.12 


1.29 


0.93 




-0.013 


-0.021 


1.68 


0.98 


6_ 


0.011 


0.016 


1.50 


0.91 



TABLE III: Form factors of B -> £>** transitions calculated in the ISGW2 model 4: 





F(0) 




a 


b 


F 


F(0) 


F(gLJ 


a 


b 




0.18 


0.24 


0.28 


0.25 




0.18 


-0.008 






A BD^ 


-0.16 


-0.21 


0.87 


0.24 


bdT 


0.18 


0.23 


0.89 


0.25 


BD 1 ' 2 

V l 1 


-0.19 


0.006 






BD 1 ' 2 
V 1 
v 2 

bd\ /2 


-0.18 


-0.24 


0.87 


0.24 


A BD " /2 


0.16 


0.19 


0.46 


0.065 


0.43 


0.51 


0.54 


0.074 


v bdT 


0.40 


0.32 


-0.60 


1.15 


bdI 12 

v 2 


-0.12 


-0.19 


1.45 


0.83 


h 


0.011 


0.014 


0.86 


0.23 


k 


0.60 


0.68 


0.40 


0.68 




-0.010 


-0.013 


0.86 


0.23 


6_ 


0.010 


0.013 


0.86 


0.23 



general, the form factors evaluated in the original version of the ISGW model are reliable only at 
q 2 = q^, the maximum momentum transfer, because the form-factor q 2 dependence is proportional 
to exp[— (q^ — q 2 )] and hence the form factor decreases exponentially as a function of (q^ — q 2 ). 
This has been improved in the ISGW2 model in which the form factor has a more realistic 
behavior at large (q^ — q 2 ) which is expressed in terms of a certain polynomial term. 

In the infinite quark mass limit, all the heavy-to-heavy mesonic decay form factors are reduced 
to three universal Isgur-Wise (IW) functions, £(uj) for s-wave to s-wave and ti/ 2 (uj) as well as 
T 3/2( UJ ) f° r s-wave to p-wave transitions, first introduced in [37|. Specifically, the B — > D$ and 

1 /2 3 /2 

B — > form factors are related to T\/2\w), while B — > and B — ► D 2 transition form factors 
are related to t^/2(uj) |43j: 

(D(v')\V,\B(v)) = eH(« + t/) M , 
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(D*(v',e)\V, 
(B*(v',e)\A„ 
(D* (v')\A, 

(Dl /2 (v',e)\A, 
(D*(v',e)\V, 



B(v 
B{v 
B(v 
B{v 
B{v 
B{v 

B{v 

B{v 
B(v 



:i+u;)e*-(e*-v)v' 



i2r 1/2 (uj)(v - v')^ 
-t2r 1/2 (cj)[(l-w)e* + ■ v)v'^ 
-2T 1/2 (Lo)e^ af3 e* u v' a vP, (2.21) 
i-±= r 3/2 M{(l - oj 2 )s; - (e* ■ v)[3v, + (2 - w )^]}, 

-^/| T 3/2(^)(1 + uje^ape V V 1 , 

\/3 r 3 /2 (u^e^a/je*^ 7 ^?/ ?/ , 
-iV3r 3/2 (o;){(l + W ) £ ;X - e^}, 



where cv = v ■ v' . For completeness, we have also included B —* D,D* transitions in terms of the 
Isgur-Wise function £(u>). Using the vertex functions in the heavy quark limit given by Eqs. (4.23) 
and (4.24) of [3| in conjunction with HQET, we have derived the IW functions in the light-front 
model. Their numerical expressions are given by 



= 1 - 1.22(w - 1) + 0.85(w - l) 2 , 
t 1/2 (lo) = 0.31 (l - 1.18(w - 1) + 0.87(u - l) 2 

r 3/2 (u/) = 0.61 (l - 1.73(w - 1) + 1.46(w - l) 2 



(2.22) 



They are in good agreement with the lattice results Ti/ 2 (1) = 0.38 ± 0.05 and t 3 / 2 (1) = 0.53 ± 0.08 
Q- 8 

It is easily seen from Eq. (|2.21|) that the B — > D** matrix elements of weak currents vanish at 
the zero recoil point u = 1 owing to the orthogonality of the wave functions of B and D**. From 



Eqs. (|2~ToT) and (|2~2T|) it is clear that the B 
limit are related to Ti/ 2 (ct>) by 



T1/2H 



an d -B — ► D^ 2 form factors in the heavy quark 
^m B m D * 



f bd Hq 2 ) 



m B - m D * 



m B - m D * 



1 



(m B ~m D *y 2 



m B m D i/2 



A BDl" 



m B - m 1/2 



m B m D i/2 



m B m D i/ 2 



T BD 



1/2 



m B - m 1/2 
^1 



V 1 (<r) 



1/2 



^ 1 (^) 



m B m D i/2 



m B - m D i/ 2 



1/2 



O (q 2 ) 



( m B-m fl i/2) 2 



(2.23) 



BD l/2 

Hence, in the heavy quark limit, we have V 1 (q ) 



vf^M? 2 ) = -A BD i 1/2 (g 2 ). Likewise, 



Comparison with other model calculations of Ti/%(1) and T3/2(l) is summarized in Table XV of |4j|. For 
example, based on QCD sum rules, t 3 / 2 (1) = 0.74 ± 0.15 is obtained [filj . 
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the B -> L>? /2 and B 



r 3/2 (' 



I?2 form factors are related to t 3/ / 2 (w) via 
~~ 2 l 3/2 (q 2 ) _ 



\ 



m B m 3/2 or 



2m B 

m 3/2 



1 

3/2/2 



\ 



2m B 

m 3/2 
1 1 



(cf(^) + c 3 _/ 2 (^) 



lm B m D * 



I m B k(q 2 ) 
3m,D* 1 + u! 

m B m D * 



2 

2^/2 



% 2 ) 



l + w 



^m B m D 3/2 q 3 /2{q 2 



with 



4/2^) 



m1)/(2m B m D * 



b-(Q 2 ) / 
■ 2 )+b-{q 2 ) 



(2.24) 



and 



RD 3/2 

■(ms-m^)^ 1 (g 2 ) 



V 



me — m D 3/2 



c 3/2 (. 2 



93/2 (? 2 



-2m D 3/2 



(q 2 



m B — m 3/2 



(2.25) 



One can check that the signs of various form factors in Table [H] are in agreement with the heavy 
quark limit behavior of B —* D** transitions, Eqs. IJ2.2HJ) and (j2.24j) . 

It turns out that, among the fourteen B — ► .D** form factors, while the covariant light-front 



l i t , ■ r 4 R n 1 /2(3/2) fin 1 

model predictions tor A i , \/q 



1/2 



1/2 



,h,b+,b- are in good agreement with those in 



.BR 



the heavy quark limit, the predictions for i^i^ ,^ 1 and k at zero recoil show a large 

deviation from the HQS expectation. Indeed, Eqs. (|2.23l) and (|2.24j) indicate that except for 
F l , these form factors should approach to zero when q 2 reaches its maximum value, a feature 
not borne out in the covariant light-front calculations for finite quark masses. This may signal 
that Aqcd/^-q corrections are particularly important in this case. Phenomenologically, it is thus 
dangerous to determine all the form factors directly from the IW functions and HQS relations since 
1 /itlq corrections may play an essential role for some of them and the choice of the [3 parameters 9 
for s-wave and p-wave wave functions will affect the IW functions. 



1/2(3/2) 



III. B D**ir DECAYS 

Given the decay constants and form factors discussed in Sec. II, we are ready to study the B 
decays into p-wave charmed mesons. In this section we will focus on B — > D**ir decays. The exper- 
imental results for the product of the branching ratios B(B — > D**tt) and B(D** — > two particles or 
three particles) are summarized in Table ITVl 



9 is a wave function parameter which governs the behavior of the phenomenological meson wave functions, 
4> oc exp(— |p 2 |/2/3 2 ). It is expected to be of order Aqcd- 
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TABLE IV: Experimental branching ratio products (in units of 10 4 ) of B decays to D**tt, where 
D%,D[,Di,D% stand for the charmed mesons D£(2400), D[(2430), D 1 {242{)) and £>|(2460), re- 
spectively. The Cabibbo-suppressed mode B — ► D 2 + K~ is also included. 



Mode 






BaBar [53, 54] 


Belle [4, 52] 


Average 


B(B~ - 


-> D* °ir-)B(D* ° 






6.1 ±0.6 ±0.9 ± 1.6 


6.1 ±1.9 


B(B- - 


-* Dfn-)B{Df - 


-> D*+vr-) 




5.0 ± 0.4 ± 1.0 ±0.4 


5.0 ± 1.1 


B(B- - 


-* D?tt-)B(DS -» 






< 0.06 


< 0.06 


B(B- - 


-> D®tv~)B(D® -» 


D°7r-vr+) 




1.85 ±0.29 ±0.35^46 


1 ok+0.45 

J - 8b -0.65 


B(B- - 


-> Z^vr- -> 




5.9 ± 0.3 ± 1.1 


6.8 ±0.7 ±1.3 ±0.3 


6.2 ±0.9 


B(B- - 


-> D^ 7T~)B(-D^° 






< 0.22 


< 0.22 


B(B- - 


-> D| 7T-)H(-D|° 




1.8 ±0.3 ±0.5 


1.8 ±0.3 ±0.3 ±0.2 


1.8 ±0.4 


B(B~ - 


-> -D2 71" )^(^2 


7-1-1- — \ 

-> D + ) 


2.9 ± 0.2 ± 0.5 


3.4 ± 0.3 ± 0.6 ± 0.4 


3.1 ± 0.4 


B{W - 


+ ^ + 7T-)^(^ + 


— ► D°7r~) 




0.60 ± 0.13 ± 0.15 ± 0.22 < 1.20 


< 1.20 


B{W - 




-» L>* tt + ) 




„ in on 

0.14 ± 0.13 ± 0.12±g;?g < 0.70 


< 0.70 


B(lf - 




» D* + vr"7r + ) 




< 0.33 


< 0.33 


B(B° - 




* D + 7r~7r + ) 




89 ± 15 ± 17 + S'22 




B(lf - 


- Dt^)B{D+ - 


♦ J D*°7T+) 




3.68 ± 0.60±8^g±gS 


7+1.1 
— 0.8 




+ D* 2 + 7T-)B(D*+ 






< 0.24 


< 0.24 


B(B^ - 


- D* 2 + tt-)B(D* + 


-» L>* 7T+) 




2 45 + 42+°- 35 + - 39 

It U.4Z_q 45 _ 1? 


2 4 +0 - 7 


B(jf - 


-» D* + 7r-)H(L»*+ 


-» Z?°7r+) 




3.1 ±0.3 ±0.1+^ 


q 1+0.4 
°-- L -0.3 


B(lf - 


-> D* 2 + K-)B(D* 2 ^ 


-► L>°7r+) 


0.183 ±0.051 




0.18 ±0.05 



TABLE V: Experimental branching ratios for B — > D**it decays (in units of 10 4 ), where I?** 
denotes a generic p-wave charmed meson. 



Mode Expt 


Mode Expt 


B~ -» D^ 7r- 9.2 ± 2.9 
B- Dfvr- 7.5 ± 1.7 
5- -» D?tt- 13.5i^ 
5- _> D| vr- 7.4 ± 0.8 


if ^ L><5 + 7T- 0.90 ± 0.45 < 1.8 
t^H-tt- 0.21±g"ig < 1.1 

^ £,1+^- 8.3+1-2 



To determine the absolute branching ratios for B — > D ir, we need some information on the 

1/2 

branching fractions of D**. The decay Dq undergoes an s-wave hadronic decay to Dn, while 
can decay into D* by s-wave and d-wave pion emissions but only the former is allowed in the heavy 
quark limit m c — ► 00. Hence, we shall assume that the Dq and D[ widths are saturated by Dtt and 
D*n respectively, so that 

B(D* ° - D+7T-) = I B(D[° - D*+7T-) = I (3.1) 

In the heavy quark limit where the total angular momentum j of the light quark is conserved, s- 
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wave D x — » I?7r is prohibited by heavy quark spin symmetry. Therefore, for Di(2420) we assume 
that the dominated strong decay modes are (D\ — > D*7r)d- W avc> (-D7r7r) p _ wave , (-D*7r7r) p _ wa ve- From 
Table ITV1 it is clear that among the possible strong decays of D±, the three-body mode D*tttt is 
suppressed relative to Dtttt. Moreover, the analysis of Z?i(2400) — > L>7r + 7r _ by Belle [5^ shows 
that the decay mode D\ DqTt gives the best description. Therefore, 

B(£>? - DV+tt-) « |b(£>° - D*+tt-), - D* + 7r-, D^tT) = |. (3.2) 

The tensor meson D| decays into D* or D via d-wave pion emission. Since the production of D*7rn 
in £>2 decay is very suppressed, we take 

B{Df -» D*+vr-, D+vr-) = H (3.3) 

In heavy quark effective theory, it is expected that 

T(Df -> D+TT-) _ 2 m D f p c (D* 2 -> Dtt) \ 5 _ 

r(D*° D^+vr-) 3 m D * \p c (Dl -> D*tt) J ' l ' j 

in excellent agreement with the direct measured value of 2.3 ± 0.6 ji^]. Applying Eqs. (|3.1|) - (|3.3jl . 
the absolute branching ratios of B — > D**ir are shown in Table [V] Note that in the factorization 
approach it is expected that 

B(B° -> ZJ 2 * + K-) . 2 

R =sm(/p, 3.5) 

_> D*+ir-) ° 

with 0c being the Cabibbo angle. From Table ITVl we see that this relation is well satisfied experi- 
mentally. 

We shall study B — * D**ir decays within the framework of generalized factorization in which 
the hadronic decay amplitude is expressed in terms of factorizable contributions multiplied by the 
universal (i.e. process independent) effective parameters a\ that are renormalization scale and 
scheme independent. Apart from a common factor of GpVcbV* d /y/2, the factorizable amplitudes 
for B- -» D**°tt- read 

- £>5°0 = - ai f n (ml-ml )F^(ml) + a 2 f Do (ml-ml)F B ^m 2 Do ), 
A{B - ^ d'^-) = 2 (e* ■p B ){a x UV^{ml)m D > i - ^^(m^m^ j, 

A(B- - D?tt-) = 2(e* •p B )|ai/ 7r y BDl (m2)m Dl - ^/^^(m^Jm^j, 
A(B~ -» D*°7r-) = -axUel^B \k(0 + &+«)(m| - m^) + , (3.6) 



with 



v m D ^ = Vq rn D i/2 cos + V m D s/2 sm0, 

R/~) BD^ BD^ 

Vn 1 fno 1 = -V 1 m 1/2 sin + Vq 1 m 3/2 cos9, 

fD'^D^ = f D i/2m D i /2 cos0 + / D 3/ 2 m D 3/ 2 sin0, 

fD 1 m Dl = - f D i /2 m D i/ 2 sin6 + f D 3 /2 m D 3 /2 cos 9. (3.7) 
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The decay amplitudes for B° — ► D** + tt~ can be obtained from A(B~ — > D**°7r~) by setting 02 = 0. 
10 Note that except B — » D 2 °tt~ all other decay modes receive contributions from color-suppressed 
internal VF-emission. In the heavy quark limit, 

-> Do°vr _ ) = -2a 1 f 7Ty /m B m Do (m B + m Do )(uj - \)t 1/2 {lo) 
+2a 2 f Do Ft(m 2 Do )(m 2 B -ml), 

A(B~ -> DfO = 2(e* • a x f^^m^m^{m B - m^/^T^u) 



-02 f D i/2Ff n (m 2 Dl/2 )m B m D i/2j, 

A(B~ -> -D?vr _ ) = \/2ai(e* • v)f n ^/m B m Dl (m B - m D 3 /2 )(w + l)r 3/2 (u;), 
A(B~ ^ Z^V") = -V^ai/ T £*,//^/mp^(m B + mfl 2 )r3 /2 (w), (3. 

where u = (rog + m^** — m 2 ) / (2m B m£)**). 
The decay rates are given by 11 



r(B-2?jp7r) = -^lifB -2?P7r)/( e *-p B )| a , 

r(B - DJtt) = -^V- (— Y 1^(5 - ^vr)| 2 , (3.9) 
l2irm z D2 \m D ,J 

where A{B — > -Djj 71 ") = £ uuP B Pb M(B — > -D2 71 ") an d is the cm. momentum of the pion. From 



Eqs. (pT5J) and (|33)) . we obtain 



r(s° - 


- ^ * + < 


) = T(B° - 






r(s° - 




= T(B° - 


♦ ^2 + 


- GF IVtV* l 2 n 2 f 2 ™ 3 
~ TEtT a lJ-K m B 


r(B~ - 


- 0?O 


= T(B° - 


» D+7T-), 


r(B~ -> Dfir-) 



(1 _ r )5 (1 + r) 3 

^ F1/2MI , 

:i-r) 5 (l + r) 7 2 

-I r 3/2(^)l , 



16r 3 

T(5° -> D^O. (3.10) 
in the heavy quark limit, where r = m^** /wig. It is evident from Table fTVl that the HQS relations 
T(B° -> D* 2 + i:-) = F(B~ -» Dfir-) and T(B° -» Dfir~) = T(B° -» L>2 +7r ~) are respected by 
experiment, while the relation T(B° — ► Dfir~) = T(B~~ — ► Di~it~) is not satisfied, implying the 
importance of color-suppressed contributions which vanish in the heavy quark limit. 

From Eq. (|3.6j) we see that apart from the coefficients a\ and a 2 , the color-allowed and color- 
suppressed amplitudes for B~ — ► {Dq°, .Df, D^lvr - have opposite signs, 12 in contrast to the case 
of B — > Z?7r decays. 



10 It is customary to neglect the ly-exchange contributions to 5° — > D** + n and £>° — > Z?**°7r°. 

11 Because the scalar resonances Dq and Z?i have widths of order 300 MeV, we have checked the finite width 
effects on their production in B decays and found that the conventional narrow width approximation is 
accurate enough to describe the production of broad resonances owing to the large energy released in 
hadronic two-bod y de cays of B mesons 0- 

12 We disagree with Q on the signs. If we follow Q to define (P{q)\A^\0) = f P q^ and (d\ /2 ( P , e)|A M |0) = 
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TABLE VI: The predicted branching ratios for B — > D** + tt~ decays (in units of 10 -4 ) calculated 
in the covariant light-front model and its extension to the heavy quark limit (denoted by HQS). 
The parameter a± is taken to be a\ = 0.88. Experimental results are taken from Table IVl 



Mode 


Theory 


HQS 


Expt 




3.1 


1.7 


< 1.8 


B° - D'+tt- 


0.8 


1.5 


< 1.1 


B° -> D+7T- 


10.4 


11.1 


7 6 +L7 
'•0-1.4 


B° -» D* 2 + ir- 


6.9 


10.8 


°- 6 -1.0 



A. Color-allowed B — > £>**7r decays 

We first discuss the color-allowed 5° — ► Z)** + 7r~ decays governed by the parameter a%. This 
is the place where the calculations are considered to be more robust. To proceed, we first fix 
a\ to be 0.88. In Table fVTl we show the predictions of B(B — > D** + tt~) in the covariant light- 
front model for B — ► D** form factors and its extension to the heavy quark limit with the IW 
functions given by Eq. ([2.22ft . 13 It is evident that the color-allowed modes are "normal". To 
illustrate this point, we consider the decay amplitudes in the heavy quark limit given by Eq. (|3.8[) . 
We see that the color-allowed a\ amplitudes for Dq and D[ production are suppressed relative to 
that for D\ and D\ production because of the smallness of ti/2( uj )/ t 3/2( uj )> (uj — l)/(w + 1) and 
(m b — m D^)l [m b +niDi)- Note that the first three modes in Table IVII prefer an a\ smaller than 
unity, whereas the Z?2 +7r_ channel favors an a\ close to unity 

If we treat the B — > D$ , D[ , D\ transition form factors as unknown parameters, we can determine 
them from the data. In order to satisfy the constraint B(B° — ► Dq + it~) < 1.8 x 10~ 4 , it follows 
that the B — ► Dq form factor is constrained to be 

Fo D Ho) £0.18. (3.11) 

This is smaller than the CLF prediction, F o(0) = 0.24 (cf Table ITTj) . From the measurements of 
B" ^ D[ + tt-,D+ ir we obtain 

V BD,1 (0) < 0.15, V BDl (0) = 0.39±g;g|. (3.12) 

BD f 

Taking V x (0) = 0.15 and using the experimental central value for the D[ — D\ mixing angle 
9 = 5.7° [Eq. J23J>], we find 

^o BDll/2 (0) = 0.11±°;°°, V BD " /2 (0) = 0.41 ± 0.04, (3.13) 



1 /2 

—jD 1 /2 rn D 1/2 s l i for decay constants, we find that the matrix elements (Dq\A^\B) and (D/ |V^|£?) should 
have signs opposite to that given in Eq. (28) of [2^. 

The predicted rates for Dq + tt~ and DI^-k - are somewhat different in the covariant model and its heavy 
quark limit extension. This is because the form factors F BD °(q 2 ) and k(q 2 ) do not respect the HQS 
relations (|2.23|1 and (|2.24(l satisfactorily. 
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to be compared with the CLF model predictions (see Table HH): V 1 (0) = 0.08 and V 1 (0) = 
0.47. 

The IW functions Ti^iyj) and t 3 / 2 (<jj) can be extracted from the data of B° — > D** + ir~: 
B° DI + tt- => |r 1/2 (1.36)| < 0.22, 

5° Z^+tT |r 1/2 (1.32)| < 0.19, (3.14) 

if ^ £>Vf tt~ => |r 3/2 (1.32)| = 0.30 ±0.03, 

B° - D*+n~ => |r 3/2 (1.31)| = 0.33^ . 

As stressed in passing, the HQS relation r(B -> £q + 7r~) = r(S" -» £>?7r - ) is badly broken and 
hence t 3 / 2 (u;) cannot be reliably extracted from Dfn~ production. Our predictions |43j 

r 1/2 (l) = 0.31, 7i /2 (1.32) = 0.22, r 1/2 (1.36) = 0.21, 

r 3/2 (l) = 0.61, r 3/2 (1.31) =0.37 (3.15) 

are in good agreement with the phenomenological determination (|3.14|) and the lattice calculations 
[5oj | . For comparison, the phenomenological determination of IW functions in 26] is given by 

r 1/2 (l) < 0.26, r 1/2 (1.32) < 0.20, 

t 3/2 (1) = 0.46 ±0.18, r 3/2 (1.31) = 0.31 ±0.12. (3.16) 

In short, Eqs. (|3.11j) . (|3. Ilf) - (j3.13f) and 1)3.14(1 are the main results in this subsection. 



B. Class-Ill B~ -> D**°ir- decays 

We next turn to the so-called class-Ill decays B~ — > D**°7r~ that receive both color-allowed 
and color-suppressed contributions. The experimental observation that the production of broad 
D** states in charged B decays is more than a factor of five larger than that produced in neutral B 
decays fTable lVj) is astonishing as it is naively expected to be a factor of two difference at most. For 
example, Dq°tt~ and Dq + tt~ rates are predicted to be similar in , while D'®tt is predicted 

to be even smaller than D^ + 7r~ in 24]. The constructive interference in B~ — ► {Dq°, D'i}tt~ 
decays requests that the relative sign between the real parts of a\ and a 2 be negative as noticed 
in passing. Moreover, if we take the Belle measurements of B — ► {Dq, D[}ir~ seriously, it will 
imply a color-suppressed contribution larger than the color-allowed one, even though the former is 
1/rriB suppressed in the heavy quark limit ! Since the color-allowed modes have been shown to be 
"normal" before, anything unusual must arise from the color-suppressed contributions. 

Before proceeding, we need to specify the a 2 parameter. Recall that lad = 0.88 ± 0.06, 
|a 2 | = 0.47 ± 0.06 and a 2 /a 1 = (0.53 ± 0.06) exp(i59°) are obtained in |3l| by a fit to the 
data of B —> Dir. Owing to the missing B -» D**°ir° decays, one can only determine 
|1 — xa 2 /ai| from the measurements of B~ — > D**°-k~ and B° — > Z)** + 7r - , where x(Dq-k) = 
/D D (m| — m^)FQ^{m 2 D ^)/[fT T {rn? B — mjj )FQ D ° '(m^)], for example. That is, a determination of 
the relative strong phase between a\ and a 2 has to await the measurement of the neutral mode 



16 



B° — ► L>** + 7r . For the present purpose, we shall choose a\ = 0.88 and a<i = —0.47 without 
considering their relative strong phase. Using this set of effective Wilson coefficients, we obtain 

B(B~ -» Dfir~) = 7.6 x 10" 4 , B(B° -> D 2 * + = 6.9 x 10" 4 , (3.17) 

for B — > D^tt^ decays, which are in agreement with experiment. 

In order to accommodate the experimental observation B(B~ — ► Dq°tt~) > 5B(B° — ► Dq + 7t~), 
the decay constant of -Dq cannot be too small. A fit to the B~ — ► Dq°it~ rate yields 

/d 5 = 148; 4 °MeV, (3.18) 

where F ( ^ Z)o (0) has been set to 0.18 [see Eq. (|3.18[l ]. This value of the Dq decay constant is larger 
than our CFL prediction f D * = 86 MeV [cf Eq. (|2.12j) ]. Putting the form factors (|3.12|) back to 

Eq. 1)3 . 6 j) (with V x (0) = 0.15) and fitting to the measured rates of B~~ — > D'iTt~ , D^-k^ give rise 
to 

f D , = 151±§5 MeV, / Dl = 73t| MeV. (3.19) 

Applying the experimental central value for the D[ — D\ mixing angle 9 = 5.7°, we find 

f D t /2 = 1431™ MeV, / D 3 /2 = 88l^ 4 MeV. (3.20) 

Eqs. (|3.18|) - (|3.20j) are t ne main results in this subsection. 

It should be stressed that if | a-2 1 is chosen to be smaller, say 02 = —0.30, then the decay constants 
/d , fjjr and all have to be scaled up by a factor of 0.47/0.30. This will lead to a decay constant 
of Dq larger than that of the D meson. This is in contradiction to the observtion that fjj* should 
be smaller than fjj (see Eq. 1)2.14(0 . Hence, |o2/oi| is preferred to be larger than naive expectation. 

Now we can get into more detail about the relative strength of color-allowed and color- 
suppressed amplitudes in B~~ — > {Dq°,Di}it~ decays. Since /d ,/d' Si /td (wd ) = 0.30 3> 
Fo Do (ml) < 0.18 and F 1 l (m D *) = 0.37 > V '{ml) < 0.15 it is evident that the color- 
suppressed amplitude is slightly larger than the color-allowed one for | (12/01 1 = 0.53. In con- 
trast, internal JU-emission is suppressed for D\ and D\ productions. Under the factorization 
approximation, the color-suppression amplitude is prohibited in B — > D^tt. It also vanishes in 
B — > D\tt in the heavy quark limit as 6 — > and / n 3/2 — > 0. Therefore, it is expected that 

T(B° -» £>2* +7 = r(B - -» Df 71 "") in general and T(B° -> Df7r - ) = r(B~ -> D?7i-) in the 
heavy quark limit. 

Comparing Eq. (|3.20|) with Table [H] and Eq. (|2.12j) , it is clear that form factors and decay 
constants extracted from the data are consistent with the CLF model calculations except for the 
decay constant / n 3/2 which needs to be positive. This can be seen from Eq. (|3.6|) that in order to 

account for B{B~ -> DJtt - ) > B(E -> D^tt - ), it is necessary to have a constructive interference 

between color-allowed and color-suppressed TV-emission amplitudes. This in turn implies a positive 

decay constant for D x ' , i.e. / n 3/2 > (This is most obvious in the heavy quark limit where — ► 0.) 

_ 1 

It is not clear to us why our light-front model prediction with a negative / n 3/2 does not work. Of 
course, it is possible that 02 is process dependent and for some reason 02/01 becomes positive for 
B~ -» D%ir~. Then / 3/2 will be negative. 
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It is worth mentioning that the ratio 



D _ B(j?-->£>3(2460)°7r-) 

i< - J B(fi-->D 1 (2420)07r-) 1 j 

is measured to be 0.54 ± 0.18 by Belle In soft-collinear effective theory, the equality of 

branching ratios and strong phases 

B(B^D* 2 M) = D * M = DiM 

B(B — > D\M) Y v ' 

holds in the heavy quark limit for both color-allowed and color suppressed modes, where M is a 
light meson - The early prediction by Neubert 0] yields a value of 0.35. 

—0 1 "5° , n'O, 



C. Color-suppressed B — * .D** u 7r u and B —* D±lu decays 

The factorizable B° — > D**°ir° amplitudes are given by 

A(B° -> D *V) = -^f Do (m 2 B - ml)F*«{ml ) : 
A(B°^D?ir ) = V2a 2 f D[ F^(m 2 D ,)m D[ (e* -p B ), 

A(B° - D^tt ) = V2a 2 f Dl F^(m 2 Di )m Dl (e* - VB ). (3.23) 

They satisfy the isospin triangle relation 

A(B° -> D^+tt-) = V2A(B° -► D**°7r°) + -> D**°7r). (3.24) 

Assuming no relative strong phases between Z)** + 7r - , D**°ir and _D**°7r° amplitudes (i.e. these 
three amplitudes are parallel or antiparallel to each other) and using the experimental data from 
Table llVl then the above isospin relation leads to 

B(B° -» D* °tt°) = (1.9 ± 1.0) x 10" 4 , B(B° -» L^V) = (2.4 ± 0.9) x 10~ 4 , 

B(B° -» D?7r°) = (3.0 ± 2.9) x 10~ 5 . (3.25) 

Since in reality there should be some relative strong phases between the aforementioned three 
amplitudes, the above predictions for color-suppressed modes can be considered as the lower bounds 
and should be robust. At the 90% confidence level, we have 

B(B° -> D* \ ) > 0.6 x 10" 4 , B(B° -► Dfir ) > 1.1 x 1(T 4 . (3.26) 

Using a2 = —0.47 and the decay constants given in (|3.18|) and (|3.19|) . a direct calculation of the 
branching ratios of the neutral modes yields 

B(B° -» D*\°) = (1.4±g;?) x 10" 4 , S(S° -> Dfvr ) = (1.4 ± 0.5) x 1(T 4 , 

£(5° -» D'jV ) = (3.2^-i) x HT 5 , (3.27) 

M24 
0. 



14 The early Belle result R — 0.77 ± 0.15 |4j does not take into account the contribution from Di(2400) 



Dir+ir-, so are the results 0.80 ± 0.07 ± 0.16 obtained by BaBar |5J and 1.8 ± 0.8 by CLEO 
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which are similar to the model-independent results (|3.25|) derived from isospin argument. It is 
important to measure these modes to see if the production of the broad D**° states in neutral B 
decays is not color and 1/rrib suppressed. Also the isospin relation Eq. Q3.24JI will enable us to 
determine the relative strong phases in B —* D**tt decays. 

Two remarks are in order: (i) Contrary to B — ► Dir decays where D°tt° rates are suppressed 
by one order of magnitude compared to D + ir~, the D**°ir rates here are comparable to D** + ir~ 
for broad D** states as the color-suppressed amplitude is larger than the color-allowed one. (ii) 
Although the B° decay into D|°7r° is prohibited under the factorization hypothesis, nevertheless 
it can be induced via final-state interactions (FSIs) and/or nonfactorizable contributions. In soft- 
collinear effective theory, this decay receives a factorizable contribution at the leading nonvanishing 
order in A/tbb p7 1. 

Although the class-Ill decays B° — > D**°ir have not been observed so far, there exists one 
neutral mode 5° — > Df lu that can be inferred from a recent measurement of 5° — ► D* + ujit~ 
decays by BaBar |58j |. There is an enhancement for D*ir masses broadly distributed around 2.5 
GeV. Assuming that the enhancement is actually due to B — > D[lo followed by D[ — > D* + ir~ 
BaBar obtained 



58] 



B(I? -> Dfuj)B{D' x -» D* + tt-) = (4.1 ± 1.2 ± 0.4 ± 1.0) x 10~ 4 . (3.28) 
It follows from Eq. (|3.1|) that 

B(B° -> Dfoo) = (6.2 ± 2.4) x 10" 4 . (3.29) 
Theoretically, the decay amplitude of ~ff° — > JDf w is given by 

^[5° Df^D' ,?>£,/ )w(e w ,j> w )] = -^e^e*"^!^ + S 2 {pb)^{pb)v + iSse^app&jfi], 

with (apart from a common factor of G^V^T/j^/v^) 

51 = a 2 (m B + m w )m D // D / Af^m^,,), 

5 2 = -2 ^ m D; / D ;AfK ; ), (3.30) 

5 3 = -2 a -l m D ,J D ,V B "{m 2 D ,). 

Then the helicity amplitudes -Hqi H+ and can be constructed as 



H, 



o 



2m D rm u 

H± = Si± rriBPcSz. (3.31) 

The decay rate reads 

T(B^D[u) = -El-(\H \ 2 + \H + \ 2 + \H_\ 2 ), (3.32) 

It is found 

B(B^ -» Dfw) = 2.6 x 10" 4 , (3.33) 

where we have assumed that B —* w form factors are the same as that for £? — ► p transitions 
which we took from 43]. The above branching ratio prediction is slightly smaller than the BaBar 
measurement. 
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TABLE VII: Experimental branching ratio products (in units of 10 4 ) of B decays to D**ir, where 
D* ,D' sl stand for the strange charmed mesons D* (2317) and 1)^(2460), respectively. 



Mode 






BaBar [59] 


Belle [60] 


Average 


B(B~ - 


-> D* s0 D°)B(D*u - 






< 6.6 


< 6.6 


B(B~ - 


-* D*qD°)B(D*q - 




10.0 db 3.0 db l-0±|;g 


9.8±H ±2.9 




B(B~ - 


-> D* s0 D*°)B(D* s q 




9.0 ± 6.0 ± 2.0±|;g 


< 9.3 


9-01^ 


B(B~ - 


-» D'^D^BiD'sJ - 






< 2.7 


< 2.7 


B(B- - 


-» D' sl D Q )B{D' sl - 






< 2.6 


< 2.6 


B(B- - 


-» D'-.D^BiD'-, - 


» ^rvr ) 


27 ± 7 ± 5+J 


11.6+jj±3.5 


14 1+ 4 ' 7 


B{B- - 


- L££>°)B(Z£ - 






< 5.8 


< 5.8 


B{B- - 




* Dsl) 


6.0 ± 2.0 ± 1.0±f;g 


5.9±{^ ± 1.8 


5.9l1i 


B{B- - 


- l£l>«°)B(I>£ - 




76 db 17 ± 18+.^ 


22.3±|;f ± 6.7 


28.0i}°; 5 7 


B(B~ - 


- l££>*°)B(I>£ - 




14.0 d= 4.0 d= 3.0±|;g 


9.8+J; 4 , ± 2.9 


11 1+ 3 - 7 


B(lf - 


+ z?: -/j+m/j: - - 






< 10.3 


< 10.3 


B(B° - 


-» z^-d+m^v - 




18.0 ±4.0 ± 3.0+1 [j 


10.3+j;g±3.1 




B(B° - 


- ^: - J D*+)e( J D: - 


- ^-°) 


15.0 ± 4.0 ± 2.0±|;g 


5.6±1;| < 13.5 




B(lf - 








< 2.7 


< 2.7 


B(lf - 


-» D£D+)£(l£ - 






< 3.3 


< 3.3 


n 

B(W - 


+ Z££>+)B(l£ - 


► £>r^°) 


28 ± 8 ± 5t^ 


16-9+le ±5.1 


19.3+J;° 


B(jf - 




> or 7) 




< 6.3 


< 6.3 


B(jf - 


+ l£d+)0(l£ - 


► ^-7) 


8.0 ± 2.0 ± i.o±f;o 


7.ll^±2.1 


7 4+2-1 
'•^-1.9 


jB(S° - 


- D'- 1 D*+)B(D'- 1 - 


- ^r^°) 


55 ± 12 ± 10±\l 


28.8+^6 ± 8.6 


35 n +11 - 7 

oo.u_ 10 4 


B(B^ - 


-> D' sl D* + )B(D' sl - 


- £77) 


23 ± 3 ± 3+j 


14.3+|?±4.3 


1 7 n+ 4 - 5 



IV. B -> DECAYS 

Under the factorization hypothesis, the decays 1? — ► D**D receive contributions only from the 
external lU-emission diagram, as the penguin contributions to them are negligible. More precisely, 
their factorizable amplitudes are given by 

A(B^D*;D) = ^j^'ax-Fa^ + ax 

+2%^(a 6 + a 8 )(Sr|s(l + 75)c|0)p|c(l-75)6|S)|, (4.1) 

where (9x^2) = 917/^(1 — 75)52- Since the tensor meson cannot be produced from the V — A 
current, the B decay into D* 2 D^ is prohibited under the factorization hypothesis. Except for 
D* 2 , the measurement of B — > D** D can be used to determine the decay constant of D** to the 
approximation that penguin contributions are neglected. 

The current measurements of B — * D**D are summarized in Table IVTT1 They are consistent with 
relation T(W -> D**~D + ) = F(B~ -» D**°D + ), but it is obvious that one needs more improved 
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TABLE VIII: The predicted branching ratios for L^~{~(2460). Branching ratios in parentheses are 
with respect to the decay mode D* + ir°, so are the experimental results. 

Mode [12] f61] BaBar [39, 59J a Belle [3, 60] CLEO [2] 

D*+tt° 56% 43% 



o 



£>+7 13%(0.24) 27%(0.63) 0.337 ± 0.036 ± 0.038 b 0.63 ± 0.15 ± 0.15 < 0.49 

0.275 ± 0.045 ± 0.020 c 0.43 ± 0.08 ± 0.04 

D*+j 12%(0.22) 24%(0.56) < 0.24 < 0.31 < 0.16 

D+ir+TT- 11%(0.20) 7%(0.16) 0.077 ± 0.013 ± 0.008 < 0.13 < 0.08 

D *sol 7% (0.13) 0.05%(1.2 x 10~ 3 ) < 0.25 < 0.58 
D*+ir 0.002% 



"The BaBar results are for the branching ratios with respect to D' sl — > D s Tr°-f arising from D*ir° and D* j. 
The BaBar data are consistent with the decay D' sl — > D s ir a j proceeding entirely through D*ir°. 
from continuum e + e~ — > cc. 
c from _B decay. 



data to test the above relation. 

Since -D* (2317) is below the DK threshold, the only allowed hadronic decay is the isospin- 
violating one, namely, D^n . Thus far, no other modes have been observed (see e.g. 39]). The 
allowed radiative mode is D* — > D*j, which is constrained to be T(D* — > D*j)/T(D* — > D s tt ) < 
0.059 by CLEO 0- Therefore, 0.94 < B[D* (2317) -> D s vr°] < 1.0. It follows from Table IVTT1 and 
(jlTTjl that 

f 58 ~ 83 MeV from B~ decays, 
aifD* = S —o (4-2) 

s0 \ 63 ~ 86 MeV from B decays. 

Hence, our prediction Jd* q = 71 MeV is in good agreement with experiment. 

Various model predictions on the branching ratios of -0^(2460) are shown in Table IVlllI In 
general, the models 61 ] and differ in the predictions of the branching fractions for the hadronic 



and radiative decays of D' sl which in turn will lead to different absolute branching ratios for B — > 
D' sl D. Very recentlyBaBar has been able to measure B{B — > D' sl D) without any assumption on 
the decays of D' sl [62j. What BaBar have measured are the decays B —> D meas Dx with -D mea s 
being a fully reconstructed charmed meson. The mass and momentum of the Dx are then inferred 
from the kinematics of the two-body B decay. By selecting final states with Dx = -Csi(2460)~, 
BaBar measurements of B(B — ► D' sl D^) are summarized in Table llXl Combining with the BaBar 
results for branching ratio products of B decays to D' sl ir fTable lyTI)) . the absolute branching ratios 
of D' sl (24Q0) are determined to be [62 1 

B(D' sl — >• D*7T°) = 0.56 ± 0.13 ± 0.09, B(D' sl -»■ D a i) = 0.16 ± 0.04 ± 0.03 . (4.3) 

It is clear from Table IVTTTI that model is more preferred. It is interesting to note that the sum 
B(D' sl — > D*ir° + Dg'y) 0.70 turns out to be the same in both models and is consistent with the 
measurement 0.72 ± 0.17 derived from Eq. (|4.3[) . 
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TABLE IX: Branching ratios (in units of 10 4 ) of B decays to D' sl D^ measured by BaBar |& 



Mode Br 


Mode Br 


B~ -> D'^D 43 ± 16 ± 13 
5° D£ D+ 26 ± 15 ± 7 


fl - -» D^Z>*° 112 ± 26 ± 20 
S° -► D£D* + 88 ± 20 ± 14 



The decay constant of is then found to be 

fl88 + t?MeV from B~ decays, , ld . 
oi/n' = < lit _ (4.4) 

\l52j^MeV from S decays. 

For comparison, the decay constants 

aifD* > 74 ± 11 MeV, oi/d/ > 166 ± 20 MeV (4.5) 

are obtained in 63]. Comparing (14.41) with (14.21) we see that fry is about two times of fn* , 
whereas fjj* and /q/ are very similar in size [cf Eqs. 1)4. 4j) and 1)3. 19|) ]. It is not clear to us why 
the HQS relation for the decay constants is satisfied for the sector {Dq,D[} but badly broken for 
{D* s0 ,D' sl }. 



The decays B — > D** D have also been studied in 271 ] . However, the branching ratios B{B - 
D* s0 D) = (3.0 ~ 3.8) x 10" 3 , B(B -► D' S JD) = (7.2 ~ 8.9) x 10~ 3 and B{B -» = (2.5 

2.9) x 10 -2 predicted in 27] are too large compared to experiment (see Tables rVTTl and ITX)) . 

For decays B~ — > D**~D*°, their factorizable amplitudes are given by 



A(B~ ^D* s0 D*°) = 2(e*- PB )a 1 / D .m D *< D *(m^ o ), (4.6) 
A[B~ -> D' s -[(e Dgl ,p Dal )D*°(e D *,p D *)] = e^e^Sig^ + S 2 {pb)^{pb)v + i^e^^p^^.], 



with 



5 X = ai(m B + m D *)m D i f D , A 1 

si si 

5 2 = -2 £ m^/^ifKO, 

5 3 = -2 ^ m D ,/ D ,F BI3 *(m D , ). (4.7) 

We consider the ratios r(B -> D**-D*)/T(B -> for DJ* = L>* (2317) and D^(2460) 

which have the advantage of being ai fp ** i ndependent. The results are exhibited in Table IXl where 
uses of the B — > D* form factors from 
consistent with experiment. 



43j have been made. It is evident that the predictions are 



V. B — > D**K DECAYS 

The experimental branching ratio products of B° decays to D* s * + K~ and D**~Tt + are summa- 
rized in Table KH The decay B -» D** + R- can only proceed through the VF-exchange diagram, 
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TABLE X: The predicted ratios T(B -» D**D*)/T(B -> D**L>) for Z>** = L>* and D' sl . Experi- 
mental results are taken from Tables IVTT1 and ILXl 



Mode 


Theory 


Expt 


D* s0 D*°/D* s0 D° 


0.49 


0.91 ± 0.73 


D* s0 D*+/D: D+ 


0.49 


0.59 ± 0.26 


D'^/D^D 


3.6 


3.4 ±2.4 


D' s - 1 D*+/D' sl D+ 


3.6 


2.6 ±1.5 



TABLE XI: Experimental branching ratio products (in units of 10 -4 ) of B decays to D** + K~ 
and D**~ir + where D* ,D' sl stand for the strange charmed mesons D*q(2317) and D' sl (2460), 
respectively. 



Mode 






BaBar [64] 


Belle [66] 


Average 


B(B^ - 


- DfK- 


) 


0.25 ± 0.04 ± 0.04 


0.24l°;i° ±0.7 


0.25 ± 0.06 


B(lf - 




) 


0.13 ±0.03 ±0.02 


0.46t[};}i ±0.13 


0.14 ±0.04 


b(b" - 




-)B(d*+ 




0.44 ±0.08 ±0.06 ±0.11 


0.44 ±0.15 


B0 3 - 


- d'£k- 


-)B(D'+ - 


-* Dfl) 


< 0.086 


< 0.086 


B(Tf - 




-)B(D£ - 




< 0.25 


< 0.25 


B(lf - 




)*ipk - 


* Djl) 


< 0.04 


< 0.04 



whereas B — ► D** ir + does receive an external Vl^-emission contribution which is however CKM 
suppressed. From Table 1X71 it is obvious that F(B° -> D*£K~) > T(B° -> D+K~). These two 
decays can only proceed via a short-distance VF-exchange process or through the long-distance 
final-state rescattering processes B -> D+ir- -> D+K- and B -> £>5 + 7r- ~» ^so"^~- ( In 
fact, the rescattering process has the same topology as VF-exchange.) Since B(B° — ► D + 7r~) « 
2.8 x 10~ 3 » 23(1?° — > DfK~), it is thus expected that the decay -B° — ► DfK~ is dominated 
by the long-distance rescattering process. As B(B° — ► Dq + 7t~) < 1.8 x 10 -4 30], we will naively 
conclude that T{B° -> )/r(B° -> D+if - ) « r(B° -> J D* + vr-)/r(S° -> D+vr-) < 0.06, 

where we have assumed that the rescattering effects of D + ir~ — * DfK~ and Dq + it~ — ► D*q K~ are 
similar. This is obviously in contradiction to experiment. Nevertheless, if -D* (2317) + is a bound 
state of csdd 0|, then a tree diagram will contribute and this may allow us to understand why 
T(B° ->_D^K-) > T(B° - D+K-). 
For B — > D**~n + decays, we obtain 

23(5° -► £>*^7r + ) = (2.6±g;|) x 10- 6 , B(B° -► L£tt + ) = (1.5 ± 0.4) x 10~ 5 . (5.1) 

They are consistent with the experimental limits: B(B° — ► D*^"7r + ) < 2.5 x 10~ 5 and B(B° — > 
L'^Y 7r+ ) < 2.7 x 10~ 5 , where use has been made of Eq. (|4.3I) . 
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VI. CONCLUSIONS 



We have studied the production of even-parity charmed mesons in hadronic B decays, namely, 
the Cabibbo-allowed decays B -> D**tt and D**D(*\ The main conclusions are the following: 

• We have shown that the measured color-allowed decays B — > D** + ir~ are consistent with 
the theoretical expectation. However, the experimental observation of B~ — ► D**°ir~ for the 
broad D** states is rather astonishing as it requires a color-suppressed decay amplitude be 
larger than the color-allowed one, although the former is l/m& suppressed in the heavy quark 
limit. 

• In order to accommodate the data of B — ► D**ir, it is found that the real part of 02/01 with 
a sign opposite to that in B — > Z?7r decays, where ai and 02 are the effective parameters for 
color-allowed and color-suppressed decay amplitudes, respectively. This indicates that 02 is 
process dependent and the nonfactorizable contributions to color-suppressed amplitudes are 
not universal; that is, they are process dependent. 

• The decay constants and form factors for D** and the Isgur-Wise functions 7i/ 2 (u;) and 
T 3/2( UJ ) are extracted from the B — > D**ir data. The Isgur-Wise functions calculated in the 
covariant light-front quark model are in good agreement with experiment. 

• The color-suppressed modes ~B° — ► D**°7r° for broad D** states and 5° — > D / 1 °(2430)a; are 
predicted to have branching ratios of order 10~ 4 . Robust lower bounds on B° — ► {.Dq, -Df}7r° 
can be derived from the isospin argument. 

• The decay constants of D* (2317) and D' sl (2460) are inferred from the measurements of B — > 
L>**-L> to be 58 ~ 86 MeV and 130 ~ 200 MeV, respectively. The large disparity between 
fD s0 and fjj' is surprising and unexpected. The ratios T(B — ► D*q D*)/T(B — > D*qD) 
smaller than unity and r(.B — > D~ 1 D*)/T(B — > Dj-yD) larger than unity are confirmed by 
experiment. 

• The observation that the production of D*qK~ is larger than DfK~ may imply a four-quark 
structure for the scalar charmed meson D* . 

There are two main puzzles concerning B — > D**7r decays, namely, why the color-suppressed 
amplitude dominates in the broad D** production in charged B decays and why is the sign of 
02/01 different from that in B —> Dir decays ? Thus far, we have analyzed the B — > L>**7t 
and D**~D decays within the phenomenological framework of generalized factorization. In the 
QCD factorization approach, the parameter 02 is not calculable owing to the presence of infrared 
divergence caused by the gluon exchange between the emitted D** meson and the (Bit) system. 
That is, the nonfactorizable contribution to 02 is dominated by nonperturbative effects. 

In soft-collinear effective theory (SCET), the color-suppressed 5° — ► D**°tt° decay is proved to 
be factorizable. More precisely, its amplitude factors into a pion light-cone wave function and a 
B — > D** soft distribution function rather than being like the naive 02 factorization as shown in 
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Eq. (|3.23|) 57,|65j. However, it is difficult to make a quantitative prediction in SCET at this stage 
and it is not clear to us if SCET can provide an explanation of why the sign of 0*2/0,1 flips from 
B -> Dtt to B -» D**tt. 
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